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Abstract 

We present a set of exactly solvable Ising models, with half-odd-integer spin-S on a square-type lattice 
including a quartic interaction term in the Hamiltonian. The particular properties of the mixed lattice, 
associated with mixed half-odd-integer spin-(S,l/2) and only nearest-neighbor interaction, allow us to map 
this system either onto a purely spin-1/2 lattice or onto a purely spin-S lattice. By imposing the condition 
that the mixed half-odd-integer spin-(S,l/2) lattice must have an exact solution, we found a set of exact 
solutions that satisfy the free fermion condition of the eight vertex model. The number of solutions for 
a general half-odd-integer spin-S is given by S + 1/2. Therefore we conclude that this transformation is 
equivalent to a simple spin transformation which is independent of the coordination number. 
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Two-dimensional lattice models are one of the most interesting subjects of statistical mechanics, both 
experimentally [TJ [2] and theoretically. Several approximation methods are used to investigate these models on 
the lattice, such as mean-field theory [TJ [3], the Bethe approximation [4j, the correlated effective field theory [5], 
the renormalization group [6), series expansion methods [7), Monte Carlo methods [8] and cluster variation meth- 
ods. Following Onsager's[5] solution for the square two dimensional Ising lattice, other solutions for regular 
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two-dimensional lattices have been considered, such as triangular [10j[TT], honeycomb [12] [13], kagome[14] lattices 
and others were explored in several works and their importance in statistical physics has waked up the search 
for a set of completely solvable models. Some other exact results have been obtained with restricted parameters, 
as investigated by Mi and Yang[TC] using a non-one-to-one transformation [T5] . 

Some half-odd- integer spin-S Ising models were already discussed in the liter ature[17]. Using the method 
proposed by Wu[18j. Izmailian |19| obtained an exact solution for a spin-3/2 square lattice with only nearest- 
neighbor and two-body spin interactions. Izmailian and Ananikian|20| have also obtained an exact solution for 
a honeycomb lattice with spin-3/2. A particular case of solutions of these models were obtained by the method 
proposed by JosephpT]. where any spin-S could be decomposed in terms of spin-1/2. Another interesting 
method for mapping the spin-S lattice into a spin-1/2 lattice has been proposed by Horiguchi|22j. 

Mapping between models is an important tool for the study of exactly solvable models. The aim of this 
letter is to present the transformation of higher half-odd-intcger spin-S systems into a simple spin-1/2 system. 

To demonstrate this transformation we use a two dimensional mixed spin-(S,l/2) on a square Ising lattice 
with a quartic interaction Cb as displayed in figHJ First this mixed model can be mapped onto an exactly solvable 
two dimensional spin-1/2 Ising model with a quartic interaction C a , such as presented in the literature [Hi- 
Second it is also possible to transform a mixed spin lattice Cb onto an effective spin-S lattice when we consider 
the spin-1/2 as a decorated Ising model of the lattice [23] Cb- Then the classical Hamiltonian for a mixed 
spin-(3/2,l/2) lattice is given by 



with < i,j > meaning summation over nearest interacting neighbors on the square lattice, and the last sum- 
mation is performed over all spin-3/2 sites. The coefficient Kr is the nearest- neighbor interaction parameter 

(3) 

of the bilinear term; K r corresponds to the parameter of the non-bilinear interaction, where r runs from 1 up 
to the coordination number (in our case up to 4); D is the single ion-anisotropy parameter acting on spin-3/2; 
Si represents the spin-3/2 particle; whereas at corresponds the spin-1/2 particle, with two possible values ±1 
(we use these values conveniently instead of ±1/2 for all spin di). 

First we would like to obtain an exactly solvable spin-(3/2,l/2) lattice Cb (see fig.l); to this end, we write 
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the first term of the Hamiltonian (fl]) as (Kr Si + Kr ' Sf)<jj, then for spin-3/2 we have four possible values 
±1/2 and ±3/2, thus the term Si + K ( r 2) Sf can be written 

A« =±\K^±\K^\ S« = ±1/2, (2) 

4 2 ' = ± § * w ± , S t = ±3/2, (3) 

with r = 1, . . . ,4. 

In order to project the spin-3/2 onto spin a with only two possible values ±1, we impose the condition 
\A r ^ | = \A r 2 ^ | = A r . Therefore we are able to find the parameters and K r 3 ^ as a function of A r , by solving 
the system ([2]) and ([3]), we obtain following 

— A.f , ,/i *r — A.«p 

< or , in which r = 1, . . . , 4. (4) 

K (3) _ _2 a . K (l) _ _13a 
iy r — 3 ^' r — g r 

The associated Boltzmann weight for a mixed spin- (3/2, 1/2) lattice has a similar structure to that 
of the model discussed by Wu and Lin[21], the associated Boltzmann weights are given by W({<7 r }) = 
cxp(7i 1 / 2 ,3/2)- For simplicity we consider our calculation in units of — /3. Using the solution given by 
(j4|) the associated Boltzmann weights are simplified, which read as 

wt =W{+, +,+,+) = acosh(Ai + A 2 + A 3 + Aa), 
w 2 =W(+, -,+,-) = acosh(Ai - A 2 + A 3 - Aa), 
w 3 =W(+, -,-,+) = acosh(Aj -A 2 - A 3 + Aa), 
wa =W(+,+, -, -) = acosh(Ai + A 2 - A 3 - Aa), 

w 5 =W(+, -,+,+) = acosh(Ai - A 2 + A 3 + Aa), (5) 
w 6 =W(+, +,+,-) = acosh(A 1 + A 2 + A 3 - Aa), 
w 7 =W(+,+, -,+) = acosh(A 1 + A 2 - A 3 + A A ), 
w 8 =W(-,+, +, +) = a cosh(-Ai +A 2 +A 3 + Aa), 
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where a = 2(c 4 + e 4 ). 

The lattice Cb can be transformed into an effective spin-1/2 lattice £ a , as shown in fi.g.1. Then at each site 
of C a there can be eight lines; we are thus led to consider an eight-vertex model C a with weights given by ([5|) 
and which satisfies the free fermion condition [55], 

UJ1LU2 + = (J5UJ6 + U7UJ8. (6) 

Consider the following Hamiltonian, 

&i 4 (7) 

all tetra- 
hedron 

for the effective spin-1/2 Ising model C a , which is exactly solvable using the eight-vertex model[26]. Here, 
means summation over all pairs of sites over the tetrahedron (the boldface tetrahedron in Fig.l), and Jij are 
their corresponding interacting parameters, whereas Ji,2,3,4 is the quartic interaction parameter; N represents 
the number of spin-3/2 sites on the lattice. 

After transforming the Hamiltonian (fT]) into (|7|), we relate their parameters using the Boltzmann weights, 
thus to obtain 

Jo =ln (2a) + - In (wiw 2 w 3 W4W 5 w 6 wrws) (8) 

o 

8 V WQW7W2W3/ 

j2)3 4 h (B«) (10) 

8 \W7W4W5W2' 

JiA =z ln ( J I 11 ) 

8 \W5W2W3W8' 

j 41 =lm( mW5W3 ) (12) 
j M (13) 

8 \wew4w5W3 / 

^=^(— ) ( 14 ) 

8 \W7 W4W3W8/ 

1 fw 1 w 4 w 2 w 3 \ 

^1,2,3,4 =g In I I. (15) 

The partition function Z of the decorated (or mixed) model in the thermodynamic limit is related to the 
partition function Z$ v of the effective eight-vertex spin-1/2 model ([7]) by the expression Z = c~P J °Z8 V . An 



analytical expression for the free energy of the free fermion model is well known [25] and after some manipulation, 
was expressed in the thermodynamic limit by 

f = Jo + o / ddd< t> ln[2o + 2b cos(6) + 2c cos(0) + 2d cos(<9 - 0) + 2e cos(6> + 6)}, (16) 

16tH Jo Jo 

where 

a =\{w\ + w\ + w 3 + w 4 ), b =wiU)3 - W2W4, 

c —wiWi — W2W3, d =W3W4 — Wfivg, e ~wj,Wi — w^Wq. (17) 

The system exhibits an Ising transition at the critical points 

W\ + 11)2 + Ws + 11)4, = 2 max{wi,w 2 ,W3, W4}. (18) 

At low temperature the system exhibits ordered states such as the ferromagnetic state (for W\ > W2, 103,11)4), 
the antiferromagnetic state (for W2 > 11)1,11)3,11)4) and the metamagnetic state (for W3 or W4 > w\, 11)2)- 

Thus our goal is to transform an exactly solvable Ising model with spin-1/2 via the eight-vertex model[26 
into an equivalent Ising model on the lattice with spin-3/2. For this purpose we introduce an auxiliary lattice 
Cb, with mixed spin-1/2 and spin-3/2, the schematic transformation of which is displayed in fig[T] If we consider 
the spin-1/2 model as a decoration of the mixed spin model and transform it into an equivalent spin-3/2 Ising 
model on lattice, then we can conclude that there exist a transformation from the spin-1/2 Ising model onto 
the spin-3/2 Ising model, with a non-bilinear interaction and four-body interactions terms over the tetrahedron 
or quartic interaction. 

Using the first of solution of Eq.Q, we express the Hamiltonian of the transformed spin-3/2 lattice C c in 
terms of four constrained parameters and one arbitrary parameter D, thus obtaining 

i a ll tetra- 

hedron 

1372 ^ 000 784 g cy3 c 1 o 448 o3 c*3 o3 256 ^3 c*3 c*3 A / 1 c\\ 

where the parameters Jij were already defined by Eqs.(j8|)- (fl"5"|) , and these parameters are constrained. Whereas 
the parameters A\, A2, A3 and A4 are free, as well as the single ion anisotropy parameter D. Summations are 
performed as indicated in (|7|). 
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By factorizing the terms under the summation sign, we obtain 

H 3/2 =JoN + £ DSf + J2 (-',,^7 - iSf)(7 - 4Sf )) + 

+ E ^2,3,4^1^(7-0(7-0(7-0(7-0. (20) 

all tetra- 
hedron 

We can reduce the Hamiltonian (|20[) further for both solutions in Eq.Q, using the spin transformation 
<j( k \S), where this function only take two possible values ±1, for all values of S = {—3/2,-1/2,1/2,3/2}. 
These spin transformations read as 

|(7 -45 2 ); k = l 
*W(S) = { 3 (21) 

f(13 -AS 2 ); k = 2. 

We remark that the transformation given by (|21[) are the same as those obtained by Izmailian}19] using the 
method proposed by Wu[18] (but for another model, the spin-3/2 Ising model on a square lattice). As we can 
see, the previous transformation is independent of the coordination number. Therefore it may be used to yield a 
particular case of solution, for the exactly solvable spin-3/2 Ising model in a honeycomb lattice as was obtained 
by Izmailian and Ananikian|20j. 

We are now able to write the Hamiltonian (|2"D|) as 

<>^+E^ + E^M fe) + E JwS^S^- ( 22 ) 

i (^j) aI l tetra- 

hedron 



The second term of the Hamiltonian (|22p may be written using a new spin variable transformation, similar 
to that performed by Izmailian|19j. Note that there the transformation was performed only for spin-3/2 case 
and for higher spin that method could become a more complex task. 

We also remark that the model discussed by Izmailian 1 19] can be completely re-obtained, by using our simple 
transformation instead of Wu's|18j method. 

The Hamiltonian (UJ) can be extended to arbitrary half-odd-integer spin-(S,l/2); the general Hamiltonian 
has the following form 

n 1/2 , s = E (K^Si + K^Sf + K^S^-. + K^sf^+J^DSl (23) 

<i.j> i 



where K r is the bilinear interaction parameter, and Kr are non-bilinear interaction parameters of Sf<jj, 
with u = 2, ...,5/2 + 1 and r = 1, ...,4, thus the set of transformations for each half-odd-integer spin-S have 
5+1/2 unknown parameters to be determined. We extend eq.Q to the general case (|23|) by constructing a 
Vandcrmondc-like matrix, explicitly 



(I) 3 



5 5 3 



5 s 



1 \2S^ 



(I) 



(|) 2S 



5 2S 



A, 



A 



(3) 



( .d) \ 



- 1 / 



,(2) 

/If 



.4 



/ 



where the parameters A r u ^ must satisfy the following identities | A^p 



\A 



(2) 



(24) 



(S+-) 

\Ar 2 \ = A r , in order to 



project the spin-S onto spin-1/2. By inverting the Vandcrmonde-like matrix, we are able to obtain the solution 
of the algebraic system. In what follows, we show the solutions of these system of linear equations for some 
higher spin-S values. 

For 5 = 5/2, there are three solutions: 



1067 c 11 c3 , 1 c5. u _ i 
~480~ J 12 ° + To J > K — 1 



a«(5) = { 



529 a _ 5 o3 
240 ° 6 



1 c5. 

15 D ' 



1183 



<? _ 23 o3 
480 J 12 ° 



7_ o5. 

30 ° ' 



k = 2 
k = 3. 



(25) 



For 5 = 7/2, there are four solutions: 



30251 q _ 301 c3 
13440 ° 288 



61 o5 1_ 

360 126 



5 7 ; 



cr«(5) = <! 



60577 
26880 



14887 
6720 



,5' 



5 



3047 c3 , 127 c5 
2880 ' 720 



11 

1260 



5 T 



637 c3 I VJ_ o5 1_ o7. 

720 180 315 ' 



68123 



c _ 1289 c3 I 293 c5 5_ o7. 

26880' ° 576 ° ' 720 J 252 ° ' 



(26) 
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For S = 9/2, we obtain five solutions: 



5851067 
2580480 


s- 


46573 c3 1 
41472 ° ' 


7501 
34560 ' 


2924921 
1290240 


s- 


813413 c3 
725760 ° 


, 3727 
' 17280 


5868067 
2580480 


s- 


334907 r-3 
290304 ° 


, 8117 
' 34560 


5725183 
2580480 


s- 


37337 c3 1 
41472 


3593 
34560 ' 


6651283 
2580480 


s- 


506017 c3 
207360 J 


I 18341 
' 34560 



c5 _ 97 c7 j 1_ o9. u _ i 

" 6048 " ' 2592 u ' »v ± 

c5 _ 239 c7 , 17 c-9. ju _ n 

^ 15190 ^ "l" 45360 U ' ^ 



15120 1 45360 

c5 _ 213. c7 , 43 c9. u _ o (27) 
° 6048 ° ^ 90720 a ' — J v y 

C5 _ _29_ c7 , 1 c-9. U — A 

" 6048 " ' 12960 u ' ^ 

c5 _ 1237 c7 , 13 o9. ju _ K 
° 30240 u ' 12960 u ' J ' 



Using the spin transformations given by eas.(|2T[). (|25|) . (|26|) and ([27|) . we arc able to recover the transfor- 
mation obtained by JosephpT] by considering k = 1, whereas the remaining transformation are new. 

In this report we present a transformation which maps a spin-3/2 lattice, with quartic and non- linear 
interactions terms, onto an effective spin-1/2 Ising model on lattice. First this transformation was carried out 
using an auxiliary mixed spin-(3/2,l/2) square lattice with only nearest neighbor interaction term, mapped 
onto an effective spin-1/2 or spin-3/2 lattice model, depending on which spin is considered to be the decoration 
spin. Finally, a systematic way of transformation for higher half-odd-integcr spin-S is considered, inverting a 
Vandcrmondc like matrix, to obtain a families of mapping between spin-S models and spin-1/2 models, using 
the one-to-one transformation. Therefore we conclude that there exist a spin transformation, which can be 
applied to lattice models with arbitrary coordination number, including non-exactly solvable half-odd-integer 
spin-S models. We also recovered some results previously obtained in the literature |191 [2"0] . as a particular case 
of solution. 
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Figure 1: Schematic representation of mixed spin-(S,l/2) on a square lattice (Cb), a square-type spin-1/2 Ising 
model (C a ) and a square-type spin-S Ising model (C c ) 
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